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Abstract 

Applying a sparse constraint on the beam pattern has been suggested to suppress the sidelobe of the minimum 

^~~^ . variance distortionless response (MVDR) beamformer recently. To further improve the performance, we add a 

p^ , mixed norm constraint on the beam pattern. It matches the beam pattern better and encourages dense distribution 

S . in mainlobe and sparse distribution in sidelobe. The obtained beamformer has a lower sidelobe level and deeper 

^-5 ; 

QQ , nulls for interference avoidance than the standard sparse constraint based beamformer Simulation demonstrates that 

the SINR gain is considerable for its lower sidelobe level and deeper nulling for interference, while the robustness 
^^ , against the mismatch between the steering angle and the direction of arrival (DOA) of the desired signal, caused 

O . by imperfect estimation of DOA, is maintained too. 

in ■ 

^ ' Index Terms 

m ; 

00 . 

t^^ ' robust beamforming, sidelobe suppression, mixed norm constraint, sparse constraint. 

in ■ 

o: 

'^ ■ I. Introduction 

> '■ 

• 'H , Multiple- antenna systems have received a lot of attention from both the wireless industry and academia, 

c3 I because of their strong potential in realizing high date rate wireless communications in next generation 

wireless networks. A beamformer is a versatile form of spatial filtering. It uses multiple antenna systems 

to separate signals that have overlapping frequency spectra but originate from different spatial locations. 

Beamforming has become a key technique in current and future wireless communications ^. 

The minimum variance distortionless response (MVDR) beamformer has been considered as a popular 
method for enhancing the signal from the desired direction while suppressing all signals from other 
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directions as well as the background noise [[0, but its relatively high sidelobe level would lead to significant 
performance degradation, especially with the unexpected increase in interference or background noise ^. 
To enhance the robustness in the presence of array steering vector errors, doubly constrained robust capon 
beamformer used a norm constraint on the weight vector to improve the robustness |l3l. To achieve a faster 
convergence speed and a higher steady state signal to interference plus noise ratio (SINR) Q constrains 
its weight vector to a specific conjugate symmetric form. In Q, fully complex- valued radial basis function 
(RBF) network with the fully complex- valued activation function is used in MVDR beamformer to short the 
convergence period. However, high sidelobe is another drawback, which would result in deep degradations 
in the case of unexpected interferences or an increase in noise power [2J. In order to provide sidelobe 
suppression and angle mistmatch robustness for an MVDR beamformer, a sparse constraint on the whole 
beam pattern was recently proposed in |[6l, and a robust beamformer with the sidelobe suppressed was 
obtained. 

In |l6l, the sparse constraint was added equally on both the mainlobe and the sidelobe. But the expected 
beam pattern would enjoy most of the high array gains in the mainlobe, and the array gains in the 
mainlobe are in dense distribution. The constraint should encourage dense distribution in mainlobe and 
sparse distribution in sidelobe. To further enhance the performance, a mixed norm constraint is incorporated 
to match the expected beam pattern better. Numerical evaluations show that the proposed beamformer 
achieves a lower sidelobe level, deeper nulls for interference avoidance. 

II. MVDR Beamformer 

The signal received by a uniform linear array (ULA) with M antennas can be represented by an M-by- 1 
vector, x{k), the expression of which is given by 



xik) = sikMOo) + J]/3jikMej) + n{k) (1) 

where k is the index of time, J is the number of interference sources, s(^) and /3i(k) (for 7 = 1, ... , J) 
are the amplitudes of the SOI (signal of interest) and interfering signals at time instant k, respectively, 61 
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(for / = 0, 1, ... , J) are the DOAs of the SOI and interfering signals, (pi = (Ind/X) smQi, with d being the 
distance between two adjacent antennas and A being the operating wavelength [IJ, i.e., the wavelength 
of the SOI, and n(^) is the additive white Gaussian noise (AWGN) vector at time instant k; 2l{6) is the 
steering vector in the angle 6, with its m-th element expf^C^ ~ 1)^ sin^j. 
The output of a beamformer for the time instant k is then given by 

y{k) = w"x(k) = s(k)w"a{eo) + J^/3j(k)w"a(ej) + w"n(k) (2) 

7=1 

where w is the M-by-1 complex- valued weighting vector of the beamformer. 

The MVDR beamformer is designed to minimize the total array output energy, subject to a linear 
distortionless constraint on the SOI. The weighting vector of the MVDR beamformer [HI is given by 

'Wmvdr = arg min (w^R^w) , s.t.w^a(0o) = 1 (3) 

w 

where R^ is the M-by-M covariance matrix of the received signal vector x{k), and w"a{6o) = 1 is the 
distortionless constraint applied on the SOI. 

III. The Mixed Norm Constraint Beamformer 

In order to suppress the sidelobe level of the conventional MVDR beamformer, a sparse constraint on 
the whole beam pattern was suggested in f6l . Accordingly, the weighting vector of the improved MVDR 
beamformer based on a sparse constraint (SC) is given by 

Wsc = arg min (w^R;,w + 71 ||w^A||^) , s.t. w^a(0o) = 1 (4) 

w 

where 71 is the factor that controls the tradeoff between the minimum variance constraint on the total 
array output energy and the sparse constraint on the beam pattern. The M-by-N matrix A is the steering 
matrix with a^s ( n = 1, 2, ... , N ) being the sampled angles in the [ -90°, 90°], and it covers all the N 
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Steering vectors for all possible interference with DOA in the sampling range, with oro being the DOA of 
the SOI as defined in ([T]), i.e. 



1 



expOVi) 



1 



exp {jifiN) 



exp(j(M-l)^i) ■■■ exp(j(M-l)Vw) 
2nd 



(5) 



Vn 



A 



sina,,, for n = 1, • • • , A'^ 



(6) 



\P\^lp 



and ||x||p = iYjiUif) ' is is the ^p norm of a vector x.When < p < 1, the ^p norm provides a 
measurement of sparsity for x. The smaller the value of ||x||^ is. the sparser the vector x is, It means that 
the number of trivial entries in x is larger [|6l. When p = 1, dH) is an second order cone programming 
(SOCP), and can be solved efficiently. 

The optimal weighting vector indicated by dH) can be found by an adaptive iteration algorithm [[61, [|71. 
When jc = 1, a simpler way called basis pursuit [8J, can solve dH) efficiently. 

In ©, the sparse constraint operates for all the array gains w^A in all the possible values of DOA from 
-90° to 90°, i.e. it enforces sparse distribution of the array gains in both the mainlobe and the sidelobe. 
However, the array gains are not in standard sparse distribution, but in dense distribution in the mainlobe 
and in sparse distribution in the sidelobe. Here instead of the standard sparse constraint, an mixed norm 
constraint with different norms on different lobes can be added to the MVDR beamformer to improve the 
performance. It can be formulated as 



y/MNB = argmin [w^R;,w + 72 (||w^Am|L + Hw^AsHj] 

w 

s.t.w^a(0o) = 1 



(7) 



where 



Am = 



a(0_,) ••• a(0o) 



a(0+6) 



(8) 
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a(0_9o) •■■ a(e.b-\) a(9+h+i) ••• a(0+9o) 



(9) 



Am and A5 are sub-matrices of the steering matrix A. Am is composed of 2b+l steering vectors with the 
sampled angles in the mainlobe; while A5 is constituted with the rest of the steering vectors in A. The 
product w^Am indicates array gains of the mainlobe in the beam pattern, and w^A^ indicates array gains 
of the sidelobe. The width of the first block is 2^7+1 corresponding to the mainlobe; and the other block's 
width is corresponding to the sidelobe. 72 is the weighting factor that controls the tradeoff" between the 
minimum variance constraint on the total array output energy and the mixed norm constraint on the beam 
pattern to shape the beam pattern, b is an integer representing the bounds of the mainlobe block. The 

minimization of w^A^ + w^Ac L is the mixed norm constraint. 

II 1 1 00 II '^ 1 1 1 

As the mixed norm constraint are used in (|7]), we name it as the mixed norm beamformer (MNB). 
To make a diff"erence, the beamformer dH) which was proposed in [6] is named as the standard sparse 
beamformer. Since the objective function of MNB is convex, the optimal Wmnb can be solved out by cvx 
[fm and SeDuMi [fHl. 

In the mixed norm beamformer (|7]), the term |[w^A5|| enforces the sparse distribution of the array 
gains in the sidelobe. That is to say, the number of the non-trivial array gains is much less than trivial 
gains'. And the term ||w^Am||^ lets the array gains in the mainlobe be dense. And that is to say, the 
amplitude diff"erence of the array gains in this area is little. Then the obtained beam pattern would be like 
this: most of the non-trivial elements (array gains) are in the mainlobe; and the rest trivial elements are 
in the sidelobe. Thus, with most of entries in the mainlobe being non-trivial, the angle mismatch would 
not seriously degenerate the performance; and with most of entries in the sidelobe being trivial, the eff"ect 
of interferences and background noise would decrease. 

To illustrate the why the ^1 norm encourage sparse distribution and the /'oo norm encourage dense 
distribution, a simple geometry is given in Fig. 1. In the figure the minimization of the £1 norm of a 
two dimensional vector, the £2 norm of a two dimensional vector and the /'oo norm of a two dimensional 
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vector are represented as a smallest rhombus, a smallest circular and a smallest exact square, respectively. 
The distortionless constraint can be represented as a line. The optimal solution for the minimum norm 
||w^A|| , { p = 1,2, +00) subject to the distortionless constraint w"a(6o) = 1 would be the tangent point 
where the line and the curve (rhombus, circular, square) meet. Therefore, for the £i norm minimization 
situation, two entries in the solution would has very different absolute values with high probability; and 
for the £co norm minimization situation, two entries in the solution would has considerable similar absolute 
values with high probability. 

The proposed MNB (|7]) fits the beam pattern better than the standard sparse beamformer in practice, 
and the performance of the proposed MNB would be improved. 

IV. Simulation Results 

In the simulations, a ULA with 8 half- wavelength spaced antennas is considered. The AWGN at each 
sensor is assumed spatially uncorrected. The DOA of the SOI is set to be 0°, and the DOAs of three 
interfering signals are set to be -30°, 30°, and 70°, respectively. The signal to noise ratio (SNR) is set 
to be 10 dB, and the interference to noise ratios (INRs) are assumed to be 20 dB, 20 dB, and 40 dB in 
-30°, 30°, and 70°, respectively. 100 snapshots are used for each simulation. Without loss of generality, 
p is set to be 1; Z? is set to be 23; and 71, 72 are all set to be 10. The matrix A consists of all steering 
vectors in the DOA range of [-90°, 90°] with the sampling interval of 1°. 

To quantify the performance enhancement by the mixed norm constraint. The SINR is calculated via 
the following formula: 

cr^Mb)"aieo)a"(eo)Mb) 
SINRib) = ' (10) 

Mb)"ii cr]aiej)a"(ej) + Qj Mb) 
where cr, and ctj are the variances of the SOI and j-th interference, Q is a diagonal matrix with the diagonal 
elements being the noise's variances. For a fixed b, WiBssib) can be obtained via Q. Then SINR(b) can 
be obtained via (flOl) . Simulations show that that from b = 1 to Z? = 35, The value of SINR increases 
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gradually at the beginning, and achieves the maximum SINR at bopt = 23. Then it drops afterwards. We 
can see that from b = io b = bopt, the SINR increases gradually with the increase of block-sparse 
constraint strength, and when bopt, the block-sparse constraint starts to mismatch the beam pattern, and 
the SINR drops. 

Fig. 2 shows beam patterns of the MVDR beamformer (O, the standard sparse beamformer (Hj), and the 
MNB (|7]) of 1000 Monte Carlo simulations. It is obvious that the best sidelobe suppression performance 
is achieved by the MNB ([7]). Among the three beamformers, the MNB ([7]) has the lowest array gain level 
in sidelobe area, and provides the deepest nulls in the directions of interference, i.e., -30°, 30° and 70°. 
The average received SINR by the MVDR beamformer (O, the standard sparse beamformer (HJ) and the 
MNB © are 1.2464 dB, 4.6289 dB and 5.8712 dB. 

Fig. 3 shows beam patterns of the beamformers that we have discussed, with each beamformer having 
a 4° mismatch between the steering angle and the DOA of the SOI ^. We can see that the MVDR 
beamformer has a deep notch in 4°, which is the DOA of the SOI. It can be explained by using the 
fact that the MVDR beamformer is designed to minimize the total array output energy subject to a 
distortionless constraint in the DOA of the SOI, so when the steering angle is in 4°, instead of 0°, the 
MVDR beam pattern maintains distortionless in 0° while resulting in a deep null in 4°. This observation 
shows the high sensitivity of the MVDR beamformer to steering angle mismatch. Comparing beam patterns 
of beamformers defined in (|4]) and ^, we can see that the MNB ^ further suppresses sidelobe levels 
and deepens the nulls for interference avoidance, and has almost the same robustness against mismatch. 
In the case of 4° mismatch, the average received SINR by the MVDR beamformer (|3]), standard sparse 
beamformer © and the MNB © are 0.0005 dB, 2.0163 dB and 3.2015 dB respectively. 

When b = 0, the MNB (|7]) changes into the standard sparse beamformer dH). For b i^ 0, the beam 
shaping property of the beam pattern is exploited. Our simulations also shows that the optimal block 
width b is between 21 to 27, when the directions of interferences are in [-90°, -15°) and [15°, 90°). 

Thus, our proposed beamformer provides improvements in terms of sidelobe suppression, nulling for 



JOURNAL TITLE, VOL. X, NO. X, MONTH YEAR 8 

interference avoidance, while maintaining the robustness against the DOA estimation errors, with respect 
to existing beamformers. 

V. Conclusion 

The proposed MNB shows superiority to the MVDR beamformer and the standard sparse beamformer. 
It outperforms in terms of sidelobe suppression, nulling for interference avoidance, while maintaining the 
robustness against DOA mismatch. 
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Fig. 1. The geometry demonstration of the £i norm, £2 norm and £„ norm in the place. 
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Fig. 2. Normalized beam patterns of the MVDR beamformer, standard sparse beamformer and the improved block-sparse beamformer, 
without mismatch between the steering angle and the DOA of the SOI. 
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Fig. 3. Normalized beam patterns of the MVDR beamformer, standard sparse beamformer and the improved block-sparse beamformer, 
with 4 mismatch between the steering angle and the DOA of the SOI. 



